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AN ANALYTICAL ESTIMATION OF THE EFFECT OF TRANSPIRATION 
COOLING ON THE HEAT -TRANSFER AND SKIN-FRICTION 
CHARACTERISTICS OF A COMPRESSIBLE, 

TURBULENT BOUNDARY LAYER 


By Morris W. Rubes in 


SUMMARY 


An analysis based on mixing -length theory is presented which 
indicates that surface blowing associated with transpiration cooling 
systems produces large reductions in both the heat-transfer and skin- 
friction coefficients for a turbulent boundary layer on a flat plate. 

The numerical results are restricted to the case of air blowing into 
air. The effects of blowing are indicated to be similar for high-speed, 
compressible flow to those for low-speed, incompressible flow. 


INTRODUCTION 


The frictional heating of the outer surfaces of high-speed aircraft 
has become a major problem in the design of these aircraft. Without 
thermal protection and at equilibrium conditions, the surfaces of these 
aircraft will begin reaching intolerable temperatures at Mach numbers 
even as low as 2. The designer, therefore, must provide thermal 
protection for his aircraft f s skin. This protection can be achieved in 
several ways; for example, by altering the aircraft* s shape to avoid 
sharp or pointed frontal surfaces, by providing a cooling system for the 
skin, and by providing a protective thermal insulating layer between 
the hot air in the boundary layer and the skin. A transpiration cooling 
system, one in which the coolant passes through small pores in the skin 
and into the outside boundary layer, has the advantage of providing both 
cooling of the skin and a protective insulating layer of coolant. It 
appears that these advantages make a transpiration cooling system most 
effective (ref. l). It is noted that the advantage of evaporation can 
also be incorporated into a transpiration cooling system. 

Much of the available work dealing with transpiration cooling is 
restricted to analyses dealing with the laminar boundary layer. The 
literature is quite extensive and references 2 and 3 represent examples 
of these analyses. Investigations which are more allied to the work 
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presented In this paper are represented by references 4, 5, and 6. 1 
These analyses treat the case of the turbulent boundary layer by divid- 
ing the boundary layer into two parts, a laminar sublayer and a fully 
turbulent outer region. The end results of these analyses, in effect, 
relate the coefficient of heat transfer under conditions of transpira- 
tion with two parameters; namely, (l) the velocity at the interface of 
the sublayer and outer turbulent portion and (2) a Reynolds number based 
on the distance of this interface from the surface with properties 
evaluated either at the wall or free -stream temperature. The rate of 
transpiration acts as an independent variable in the latter relation. 
Because it is known that both of these parameters are dependent on the 
local skin friction for the case of zero transpiration, it would be 
expected that their partial dependence on the local skin friction would 
continue even under the conditions of transpiration. In effect, then, 
the results of these analyses are modified Reynolds analogies relating 
heat transfer and skin friction when blowing occurs. A similar analysis 
comprises a portion of this report; however, it is modified to also 
include the frictional dissipation occurring in the high-speed air flow 
over the surfaces of aircraft. 

Several experiments have been performed concerning transpiration 
cooling in tubes or channels having turbulent boundary layers (refs. 5, 
7, and 8) . Although the results point out the advantages of transpira- 
tion cooling, the geometry of the tests and the low velocities employed 
in the boundary layer make the results of these tests too specific to 
apply in the general case of aircraft. The experiment described in 
reference 9, although also employing low air speeds, provides data on 
a flat plate, which is a fundamental aerodynamic shape. These data will 
be compared with the present analysis in a later portion of this report. 

From this brief review of literature dealing with transpiration 
cooling, it is apparent that there is no information about the influence 
of transpiration on the behavior of a compressible, turbulent boundary 
layer, such as exists on the surface of a high-speed aircraft. The 
purpose of the present report is to present an approximate analysis for 
determining the effect of transpiration on a high-speed turbulent 
boundary layer. Because of the many uncertainties inherent in the 
analysis-, it mainly has heuristic value. Experimental data will be 
required before the value of the analysis as a means of interpolating 
or extrapolating limited amounts of data can be assessed. 


1 A simultaneous investigation, essentially the same as is contained 
in portions of the present report, has been reported recently in an 
article by W. H. Dorrance and F. J. Dore, entitled "The Effect of Mass 
Transfer on the Compressible Turbulent Boundary Layer Skin Friction and 
Heat Transfer." Jour. Aero. Sci., vol. 21, no. 6, June 1954. 
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A 

SYMBOLS 

• 


a 

arbitrary constant 

A 

parameter, / 

J T w/ T i 

b 

arbitrary constant 

B 

l + Zli Ml 2 

parameter, - 1 

%/Tx 

C P 

specific heat at constant pressure 

c f 

m 

local skin-friction coefficient 

f 

m 

ratio, 7—, defined by equation (10) 

Cf/2 

F 

elliptic integral of the first kind 

g 

exponential term defined by equation (A33) and (A^b) in 
Appendix A 

G 

function defined by equation (6) 

h 

local heat- transfer coefficient defined by equation (A 78 ) 

k 

thermal conductivity 

k m> k n 

parameters of elliptic integral of first kind (see eqs. (A66) 
and (A68) ) 

K 

mixing length parameter, l = Ky 

l 

mixing length 

M 

Mach number 

Pr 

0 

(iCp 

Prandtl number, — — 
k 


e M c n 

turbulent Prandtl number, —■ -*■ 

e H 1 
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local heat-transfer rate per unit area 

local temperature recovery factor (see eq. (A80) ) 

Reynolds number based on distance along plate and free-stream 
fluid properties 

Reynolds number based on momentum thickness and free-stream 
fluid properties 

local Stanton number, 

PlUiCp 

temperature 

velocity component parallel to plate 

velocity component perpendicular to surface of plate 

coordinate parallel to axis of plate 

coordinate perpendicular to surface of plate 

parameter defined by equation (A 58 ) 

parameter defined by equation ( A 59 ) 

ratio of specific heats, 1.4-0 for air 

eddy thermal conductivity, defined by equation (All) 

eddy viscosity, defined by equation (A10) 

constant relating transpiration rate with local skin-friction 

coefficient, p F v w/ p . 1 . 1 ^ 

Cf/2 

momentum thickness 

viscosity 

mass density 

parameter defined by equation (A 6 o) , £ 
local shear stress 

parameter of elliptic integral of first kind 
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CO 

exponent in viscosity - temperature relationship 

Superscripts 

f 

instantaneous fluctuating quantity 
quantity averaged with respect to time 

+ 

dimensionless quantity (see eqs. (A 26 ) and (A27)) 

/V 

dimensionless quantity (see eq. (A19)) 

Subscripts 

a 

interface of laminar sublayer and outer turbulent region 

c 

coolant fluid initial condition 

« 

e 

u 2 

temperature recovery condition, T e = Ti + r — — 

2c p 

w 

surface 

1 

outer edge of boundary layer 

ANALYSIS 

Restrictions on Analysis 


The present analysis has several restrictions which have been 
imposed because of inherent difficulties or for purposes of simplifica 
tion. These restrictions are as follows: 

1. Only the compressible, turbulent boundary layer is considered 



(usual boundary-layer assumptions). 

2 . 

The boundary layer is on a flat plate placed parallel to the 
free-stream direction . 

• 

3 . 

The coolant fluid is the same as the boundary-layer fluid and 


enters the boundary layer at the temperature of the surface 
of the plate (no diffusion). 
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4 . The boundary layer is idealized to be composed of two distinct 

sections: a laminar sublayer next to the surface, and a 

fully turbulent outer layer. (Note: It is assumed that 

velocity, temperature, heat flux, and shear are continuous 
at the interface of the two layers.) 

5. In the fully turbulent outer layer, the transport of heat is 

proportional to the transport of momentum. 


Parts of the Analysis 


For relative simplicity the analysis is divided into two parts. 

The first part contains a determination through the use of mixing-length 
theory of the effect of transpiration on skin friction and heat transfer 
in a compressible turbulent boundary layer when the Prandtl number is 
unity. The assumption Pr = 1 greatly simplifies this portion of the 
analysis. The second part contains a determination of the effect of 
transpiration on the relation between the coefficients of skin friction 
and heat transfer when the Prandtl number is not unity. An expression 
for temperature recovery factor is also found in the second part. 

A combination of the results of the two parts (the product of 
when Pr = 1 and St/cf when Pr / l) allows the determination of the local 
heat-transfer coefficients under the conditions of transpiration when 
Prandtl number is not unity. To permit this combination, it is necessary 
to employ the premise that the local skin-friction coefficient is 
relatively independent of Prandtl number. For the case with no tran- 
spiration, this premise is good for a Prandtl number which corresponds 
to air (ref. 10) . It is not expected that small amounts of transpiration 
should alter the basic premise. 


Description of Method of Analysis 


This description of the analysis is intended to point out its 
salient features without introducing the complexities of the detailed 
mathematics. The mathematical details may be found in Appendix A. 

The basic equations that are used in this analysis are the Navier- 
Stokes equations representing conservation of momentum, the conservation 
of energy equation, and the continuity or conservation of mass equation. 
When all the quantities in these equations are represented as the sum 
of a time averaged quantity plus an instantaneous quantity, these 
equations can be made to apply to turbulent flow (ref. 11). When these 
"turbulent flow" equations are time averaged and the usual boundary- 
layer assumptions are imposed, the resulting boundary-layer equations 
are (ref. 10) : 
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— (pu) + — (pv + p'v’) = 0 


dx 


dy 


(i) 


-- SG 

pu 


Su / —7 rv 5u d /- du : r\ 

— + (pv + p’v 1 ) — = — [ [i pu T V T 1 

OX dy oy v oy / 


( 2 ) 



Barred quantities represent time averages, while primed quantities 
represent instantaneous values of fluctuating quantities. While the 
specific heat is maintained constant, the fluid properties of density, 
viscosity, and thermal conductivity are allowed to vary with temperature. 
Frictional dissipation is also included. 


Equations (l) through ( 3 ) are intractable at present, because of 
their complexity and because the behavior of the turbulent fluctuating 
components in shear flow is not known. Certain nonrigorous simplifica- 
tions are necessary to achieve a solution useful from an engineering 
viewpoint. 

The first simplification necessary in achieving an engineering 
solution to the problem is to assume in equations (l) through ( 3 ) that 
the variations of the dependent variables with respect to x are 
negligible compared to their variations with respect to y in the 
evaluation of the local velocity distribution. This allows replacing 
partial differential equations by much more easily solved ordinary 
differential equations, with y as the independent variable. The above 
assumption is based on experience with the low-speed, incompressible 
skin-friction problem without transpiration. It has been found that 
this assumption combined with a group of assumptions concerning the 
mechanism of turbulence and the structure of the boundary layer together 
with certain mathematical simplifications can be made to yield skin- 
friction coefficients which agree well with the experimental results. 


At this point it is necessary to define the Prandtl number as it 
is used in this report. In the conventional definition Prandtl number 
is the ratio of kinematic viscosity to the thermal diffusivity of a fluid 
where the transport of momentum and heat is promoted by molecular means. 
This Prandtl number is then a controlling variable for processes con- 
taining laminar flow where molecular transport phenomena occur. When 
the flow is turbulent, however, the large-scale transport of momentum 
and heat is promoted by the eddying motion, and the conventional Prandtl 
number does not govern this mechanism. In order to relate the turbulent 
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transport of momentum and heat transfer, it is necessary to define a new 
term which has the same characteristics as Prandtl number, except that 
properties defining the term are based on eddying motion. To distinguish 
this term from the conventional Prandtl number, it is called the turbulent 
Prandtl number. A turbulent Prandtl number of unity is equivalent to 
the Reynolds analogy. 

When the basic equations are reduced to ordinary differential 
equations, they are applied to the fully turbulent outer layer by 
omitting all the terms which are based on molecular transport. In the 
laminar sublayer, however, all the terms based on eddying transport are 
omitted. Thus, two sets of equations are employed, the solutions of 
each being matched at the interface of the two layers. Experience has 
shown that this two-layer model of a turbulent boundary layer is quite 
good as long as both the Prandtl number and turbulent Prandtl number are 
close to unity. 

At this stage, the analysis is divided into two parts. In the 
first part it is assumed that Prandtl number and turbulent Prandtl number 
are both equal to unity, which results in a direct relationship between 
the temperature and velocity. Because of this relationship, it is not 
necessary to solve the energy equation. The momentum equation is solved 
together with the continuity equation to yield the velocity (and tem- 
perature) distribution through the boundary layer with the local surface 
shear as a parameter. The velocity distribution is then substituted in 
the integral expression for the momentum thickness. As this integral 
is not readily solvable in closed form, the first term of a series 
solution is employed. The first term represents an asymptotic solution 
for the case of extremely low values of the skin-friction coefficient 


skin-friction coefficient and Reynolds number with the momentum thickness 
as the characteristic dimension. When the von Karman momentum-integral 
equation is integrated using the relationship c^Rg), there results an 

expression for the local skin-friction coefficient in terms of a Reynolds 
number using the length along the plate as the characteristic dimension. 
This latter integration is also not made in closed form, the first term 
of a series expansion again being used as an asymptotic relationship 
for small values of the skin-friction coefficient and the transpiration 
rate. The first part of the analysis, consequently, results in the 
determination of the effect of transpiration on the relationship 
c f (R x ) when Pr and Pr t are both unity. 

It is known that the local heat-transfer coefficient is dependent 
on Pr and Pr t . To evaluate this dependence, the second part of the 
analysis is performed to determine the relationship between the local 
Stanton number and the local skin-friction coefficient under conditions 
of transpiration when Pr and Pr t are not unity. This is achieved by 


and low rates of transpiration; that 




Thus, a relationship is obtained between the local 
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taking the ratio of members of the energy equation to corresponding 
members of the momentum equation. The ratio of the members forms an 
ordinary differential equation relating total energy and velocity. When 
the differential equation is integrated in closed form across the bound- 
ary layer, there is obtained a relationship between Stanton number and 
the local skin-friction coefficient. An ancillary result is an expres- 
sion for the local temperature recovery factor. 

A combination of the two parts of the analysis provides St(R x ) 
under conditions of transpiration for Prandtl numbers other than unity. 


Summary of Results of Analysis 


The results of part I of the analysis can be summarized as 


follows: 



= G 


w 



p W V W /PlU X 

c f/2 


or 



(5) 


where 
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when 


2 . _ 3 + q 

3 - a 


hn 


< 1 


sin 2 9 m = - 

3 + a 


or 


G = 


2K| .^) 


1/2 


> (6) 



|F[k n , <P n (u a )] - F[k n , q> n (l)]} 




+ a 


when 


3 + a 


sin 2 9 n = 


3 - u 
3 - a 


Note: 


a = 


3 = 


B - J_ B 2 + 4a 2 
2A 2 

B + <J B 2 + 4A 2 
2A 2 


a = 


(Cf/2) 

/ p W v W 
V Pi u i 


The term F(k/P) is an elliptic integral of the first kind. 


For the values of K, y^> and u a to be used in the above see "DISCUSSION" (note u a = u aA / c^/2 ) . 
The results of part II of the analysis can be summarized as follows: 



2 


(2-Pr t ) 



PW V W 

PlUl 





+ 


ff\ Pr - Pr VVw 

2 ) V P1U1 + 




-Pr) 




(8) 
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DISCUSSION 

Evaluation of K, u^, y* 


The end results of the analysis, as represented by equations (4) 
through (8), are entirely general until specific values of K, u^, and 
y+ are selected. This means that before numerical answers can be 
determined from this analysis, information other than that contained in 
the analysis must be employed. This information, that is, the values 
of K, u£, and yj, must be determined from experiment. Thus, the 
numerical results of this report can at best be as good as the experi- 
ments used in selecting K, uj, and y+. 

Low-speed case with no transpiration .- For the incompressible, zero 
transpiration case K, u£, and y+ can be evaluated in three ways. The 
first way is from velocity -distribution data, where K is a measure of 
the rate of increase of the velocity with distance from the surface, 
and ua and y+ are the conditions at the interface of the sublayer and 
outer turbulent portion. The second way is from Cf(Rg) data. In this 
evaluation the terms K, u£, and yj lose some of their physical signi- 
ficance in that they may be considered as arbitrary constants which 
absorb some of the inaccuracies introduced by the mathematical simpli- 
fications employed in evaluating the momentum thickness from the 
velocity-distribution data. The third way is from Cf(R x ) data where, 
again, K, Ua, and y+ are treated as arbitrary constants, absorbing 
still other mathematical simplifications. 

The experimental velocity -distribution data can be plotted in the 
dimensionless form u + vs. y + when local surface shear stresses have 
been obtained simultaneously. By choosing K, u^,, and y£ appropriately, 
the analysis can be made to fit portions of the boundary-layer data very 
well. If the quantities K, u£, and y£ are adjusted so that the inner 
portion of the turbulent layer is well represented, the analysis, how- 
ever, will not fit the data in the outer portion of the layer, and vice 
versa. This lack of agreement between analysis and experiment is not 
entirely surprising when it is recalled that the analysis is based on 
a primitive mixing-length theory and that several mathematical simpli- 
fications have been imposed. 

Before the experimental data of Cf(Rg) or Cf(R x ) can be used to 
evaluate the arbitrary quantities, it is necessary to relate u£ and y* 
analytically. This is necessary because the skin-friction data will 
only supply either K and u£ or K and y+ with sufficient accuracy. 

For the incompressible, no transpiration case, a comparison of the 
quantities evaluated in the three different manners is shown in the 
following table: 
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From velocity 




survey data 

From Cf(Rg) 

From Cf(R x ) 


( inner portion of 
turbulent layer) 

data 

data 

K 

0.400 

0.352 

0.392 

u a = 7a 

11.5 

12.6 

13.1 


xThese values were evaluated from figures of data summarized in reference 
12. It is apparent that the values of K and uj given by the alterna- 
tive methods do not agree , the extreme values of K differing by about 
lb percent and the extreme values of Ua also differing by about l4 per- 
cent. These differences, undoubtedly, are due to the inadequacies of the 
mixing-length theory and the mathematical approximations imposed. 

For the purposes of the present report, however, the chief concern 
is adequate representation by the analysis of Cf (R x ) . In f igure 1, 
there has been reproduced a figure from reference 12 showing available 
data of Cf(R x ). The curves given by the present analysis are represented 
by solid lines. The analytical line based on c-^R^) data differs from 
the c ; f(R x ) data by less than 7 percent. The analytical line, based on 
velocity-distribution data, is about 12-percent higher than the c^(R x ) 
data. From this it can be inferred that no matter which method is 
employed to obtain K, u+, and y+, use of these values to obtain c f (R x ) 
will result in answers adequate for engineering purposes. 

Low-speed case with transpiration .- There exist insufficient data 
to define precisely the effect of transpiration on K, u*, and y+. The 
data of reference 9 do not allow these determinations even though they 
include temperature-distribution surveys, velocity-distribution surveys, 
and measurements of local heat transfer. The temperature-distribution 
data were not used because evaluating K, u£, and y* from temperature - 
distribution surveys and local surface heat-transfer measurements 
requires making assumptions concerning the analogy between heat transfer 
and momentum transfer. Examination of the velocity -distribution data 
revealed that local skin-friction coefficients were more than 20 percent 
higher than the values in figure 2 for the case with zero transpiration. 
These high values were attributed to surface roughness effects. Since 
it is not known how roughness is affected by transpiration, it is 
believed that the velocity -distribution data may include the effects of 
both. Because heat transfer in turbulent flow seems to be affected by 
roughness to a much lesser extent than is skin friction (ref. 13 ), it 
was decided to employ the heat-transfer data in the evaluation of K, 
uj, and y+. As the data are not sufficiently precise to do this 
directly, it was necessary to employ the alternative approach of postu- 
lating values of K, u£, and y^ and then of comparing the end results 
with the data of reference 9* Physical considerations were employed in 
guiding these postulates. 
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From the physical viewpoint it would be expected that u a and y a 
should be consistent with the velocity-distribution curve within the 
sublayer, equation (A25) . In addition it would be expected that the 
Reynolds number at the outer edge of the sublayer, based on the thick- 
ness of the sublayer and equal to u£y+, should be nonincreasing under 
the disturbing influence of transpiration. Consequently, the two 
assumptions investigated are: 

(a) u a y+ = (13.1 ) 2 solved simultaneously with equation (A25) and 

(b) u a = 13.1 and y a from equation (A25) . 

The first of these represents the case of a constant sublayer Reynolds 
number, while the second corresponds to a Reynolds number which becomes 
smaller as the transpiration rate is increased. A comparison of the 
individual results based on these assumptions provides an indication of 
the influence the assumptions concerning the sublayer may have. To 
examine this point further, other numerical results were obtained using 
the condition that u£ = y£ = 13 . 1 , identical with the no transpiration 
case. It should be noted that this latter condition is inconsistent 
with the sublayer equation (eq. (A25)). 

Besides the question of the boundary conditions at the outer edge 
of the sublayer, there exists the problem of choosing a value of K, 
the mixing— length parameter, under the conditions of transpiration. Only 
well-defined experiments can answer this problem and, at present, the 
only recourse available is to assume it is unchanged by the transpiration; 
namely, K = 0.392. In addition to these assumptions, the conditions 
M = 0, Ty/Tx = 1, Pr = 0.72, and Pr t =1.00 are imposed on the analysis 
to make it conform with the test conditions of reference 9* 


Comparison with Data 


A comparison of the data and analysis is shown in figure 2 for the 
case where the transpiration rate, PyVy/piUi, is equal to zero. The 

local Stanton number is plotted as a function of Reynolds number with 
the distance along the channel wall used as the characteristic dimension. 
It appears that the data of run H-l correspond to a laminar boundary 
layer; the data of runs H-2 and H-3 correspond to a boundary layer which 
is laminar at the lower Reynolds numbers and transforms to one that is 
turbulent at the higher Reynolds numbers; and the data of run H-4 
correspond to a boundary layer which has become turbulent very close 
to the leading edge of the channel wall. When the data of run H-l are 
compared with Pohlhausen's theoretical result for laminar flow, it is 
noted that the data are generally about 20-percent higher than theory. 
Since Pohlhausen's theory was checked experimentally in the past by Fage 
and Falkner (ref. l4) , some doubt is raised concerning the accuracy of 
the data with laminar flow. The data corresponding to turbulent flow, 
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however, agree fairly well with the well-known empirical Colburn curve, 
reference 15 , or with a curve derived from the present analysis (note 
for case where Py v w/pi u i = 0, various postulates about the sublayer are 

identical) . The deviations shown by the high Reynolds number data for 
each of runs H-3 and H-4 are not considered. 

For the purpose of the present paper, the main point is considera- 
tion of the agreement of the data corresponding to turbulent flow with 
the results of the present analysis. It should be emphasized that for 
comparison purposes the data should actually be plotted against a Reynolds 
number in which the characteristic dimension is the effective length, of 
a fully turbulent boundary layer. The effective length is defined as 
the distance along the flat plate where a boundary layer, fully turbulent 
from the start, would have the same Stanton number as the actual boundary 
layer for which transition from laminar to turbulent flow has occurred 
downstream of the leading edge. The data given in reference 9 are 
insufficient to define precisely the effective start of turbulence, so 
it is necessary to resort to estimates. Because the Stanton number 
approaches infinity as Reynolds number approaches zero, any of the data 
points corresponding to turbulent flow having a finite value of Stanton 
number are downstream of the effective starting point of the turbulent 
boundary layer. Thus, if the effective starting point of a fully tur- 
bulent boundary layer is chosen as the point where the data indicate 
transition has just ended and the flow is completely turbulent, the 
effective -length Reynolds number so determined would certainly be smaller 
than the true effective Reynolds number. When this extreme correction 
is applied to the data, that is, 1.8xl0 5 is subtracted from the R x 
for run H-2, 9*10 4 is subtracted from the R x for run H-3, and the R x 
of run H-4 are left unchanged, it is found that the data cluster near 
the analytical result with a spread of about 20 percent as shown in 
figure 3* The results based on the correct effective starting length 
would undoubtedly lie between those shown in figures 2 and 3 an( i would 
be in good agreement with the proposed analytical method, considering 
the inherent scatter in the data. 

Before a comparison of the proposed analytical method and the data 
can be made under conditions of transpiration, it is necessary to discuss 
aspects of the variation of the transpiration rate along the surface of a 
plate. If it is assumed that both conduction within the porous surface 
and radiation are negligible, a heat balance on a portion of the porous 
surface results in the following equation: 


St(T e - T w ) 


Pw v w 

P1U1 


(T W 


T c ) 


( 9 ) 


Because St varies with x, it is apparent that if T e and T c are 
constant with respect to distance along the plate, the only condition 
compatible with a constant surface temperature is Pw v w/ PlUl = eSt 
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where e is a constant of proportionality. On the contrary, if 
PVfVy/piUx = constant, the surface temperature Ty will vary along the 

surface of the plate. 

In the present analysis it is quite inconvenient to treat the first 
of these cases exactly. This arises from the fact that the analysis is 
divided into two parts, that is, the local Stanton number is determined 
from a local skin-friction coefficient that was determined for the 
particular case of Pr = Pr t = 1 and a prescribed variation of 

p w v w /piui. Thus, to make PyVy/piU], = eSt would require an iteration 

process. For convenience, and because it does not produce a large 
error in the practical case, the relationship 


( 




Pr 


t> 


PW V W ^ 

PlUl’ 



( 10 ) 


will be assumed. This assumption implies that the Reynolds number 
dependence of the function f is negligible. It will be shown later 
that although this is not exactly true, the assumption only introduces 
an error of about ±4 percent over the Reynolds number range of the 
experiments in reference 9* When equation (10) is introduced into 
equation (9) there results: 


T e 



( 11 ) 


Thus letting PyVy/piUi = £(c f /2) results in essentially a constant 
surface temperature, where C is a constant of proportionality. 


For the case where the fluid properties are essentially constant, 
that is, low-speed flow with small temperature variations, it is possible 
to apply the analysis directly to the case where PyVy/piUi = constant. 

Some error will occur because of neglecting the effect of axial 
variations in the surface temperature on the convective heat transfer. 

For the case with zero transpiration, it is known that continuous 
variations of surface temperature have little effect on heat transfer 
in a turbulent boundary layer (ref. l6) . It is not expected that small 
amounts of transpiration should markedly alter this. Thus, the results 
provided by the analysis for PyVy/piUi = constant should be applicable 

for comparison with the data of reference 9* 
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Comparisons of the analytical results with the data of reference 9 
are presented in figures 4, 5> 6, and 7- Figures 4 and 9 correspond to 
constant PyVy/piUi, and figures 6 and 7 correspond to 

pyVy/piUj. = £(cf/2). In these figures no corrections for the effective 
starting length have been made for the Reynolds numbers of the data. 

In figure 4 the variation of Stanton number with length Reynolds 
number is shown for P^v^/pi u i = 0.002. The analytical results based 

on the assumed conditions of K, uj, and y+ are also included. It is 
surprising how well either postulate compatible with the sublayer 
equation compares with the data. The analysis with either of these 
postulates gives about the correct order of magnitude for the reduction 
in Stanton number due to transpiration and, in addition, indicates the 
larger variation in Stanton number with changes in Reynolds number. The 
differences in the analytical results and the data are shown much more 
clearly in figure 5 where Pw v w/ piUl = 0.006. In this figure neither 

analytical result based on either postulate agrees well with the data, 
however, both analytical results indicate a large reduction in St 
and a larger variation of St with R x , as is exhibited by the data. It 
is interesting that the two analytical results bracket the data. In both 
figures, however, the postulate where the sublayer is considered unchanged 
by transpiration leads to analytical results which exhibit the poorest 
agreement with the data. Another point of interest is a confirmation 
that variations in surface temperatures affect the data to a small 
degree. This is seen from a comparison in figure 4 of runs 17b and 20b 
with runs l6a, 17a, and 20a, where the differences between a constant 
surface temperature and a variable surface temperature are masked by 
the scatter of the data. The results of figure 5 further substantiate 
this point. 

In figures 6 and 7 the variation of Stanton number with length 
Reynolds number is shown for Py v w/Pi u i = £(c-f/2). In view of the 

previous results, only the analytical results based on postulates 
compatible with the sublayer are compared with the data. To agree with 
the experimental values of £ the values of £ for the two cases are 
1.04 and 1.15, respectively, in figure 6. The corresponding values of 
£ in figure 7 are £ = 2.68 and 3*34. Although the data scatter 
markedly, a comparison of the data and analysis is quite favorable when 
consideration is made of the probable effective starting point of the 
data of run H-10. 

From the above comparison of data of reference 9 with the results 
of the present analysis, it can be concluded that for the low speed, 
incompressible case, answers sufficiently accurate for engineering 
purposes can be provided by the analysis when the arbitrary constants 
of the analysis are set equal to: 
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K = 0.392 
u£ = 13-1 

and y+ is determined from equation (A25) . It will be noted that a 
postulate considering the sublayer Reynolds number variable is employed 
in preference to a postulate where it is considered constant. This is 
done from the viewpoint of providing a more conservative answer (i.e., 
less reduction in St with transpiration) than is exhibited by the bulk 
-of the data. 


Examples of Results of Analysis, Including Extensions 
to High-Speed, Compressible Flow 


Examples of the effect of transpiration on the local skin-friction 
coefficient, the ratio of Stanton number to the skin-friction coefficient, 
and the local recovery factor are shown in figures 8, 9 , and 10. These 
curves include cases of Mach number equal to 0, 2, and 4. The curves 
shown are subject to the following conditions: 


(a) 

K = 0.392 

- 

(b) 

u+ = 13.1 / — 
VT X 


; (c) 

y+ is determined from equation (A25) 


(a) 

00 

• 

O 

II 

3 



1 


For figure 8 the conditions for Prandtl number are: 


(e) Pr = 1.00 

(f) Pr t = 1.00 

In figures 9 and 10 the conditions for Prandtl number are: 


(g) Pr = 0.72 

(h) Pr t = 1.00 

In all the figures, the 
transpiration parameter 


parameter of the curves is the dimensionless 


i = 


(pw v w/pl u i) 

(e f /2) 


This condition of constant 


5 , 


implying a varying transpiration rate along the surface, was chosen for 
the examples from the practical viewpoint that it requires a smaller 
amount of transpiration air to maintain a certain maximum surface 
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temperature over a region -than does a constant transpiration rate 
system. Two surface-temperature conditions are presented: the first, 

is where the surface is at the total temperature and the second is where 
the surface temperature is at the free-stream static temperature. These 
surface temperatures are certainly the extreme limits of the range of 
temperatures likely to be encountered on the surface of a cooled aircraft. 

Conditions (a) through (c) are essentially the same as postulate (b) 
(page 14) defined previously, except that u& is based on fluid prop- 
erties evaluated at the surface temperature. This accounts for the 
temperature-ratio term in (b). Condition (d) represents the viscosity 
exponent corresponding approximately to conditions occurring in flight. 
Conditions (e) and (f) were chosen for purposes of mathematical simpli- 
fication while (g) corresponds to air and (h) implies an exact Reynolds 
analogy in the outer turbulent layer. 

It is apparent from figure 8 that transpiration has a very strong 
effect on lowering the value of the local skin-friction coefficient at 
all the Mach numbers considered and at both surface -temperature con- 
ditions. As an example, at M = 4 where T w = T x , the local skin-friction 
coefficient at R x = 10 7 is reduced to around 15 percent of its original 
value by a blowing rate of PyVy/pxUj. = 0.003 ( £ = 20) . 

From figure 9 it is seen that the effect of transpiration at all 
the conditions considered is to increase the ratio of Stanton number 
to the local skin-friction coefficient when Pr = 1. This increase is 
of a much smaller magnitude than the decrease in the local skin-friction 
coefficient, with the net result being that the Stanton number is also 
markedly reduced by transpiration. For example, at the point previously 
considered, the Stanton number falls to about 21 percent of its original 
value . 

In figure 10 there is shown a series of curves for the local 
recovery factor at M = 0. This figure is representative of all the 
other conditions considered, and the following discussion applies in 
general. From the figure it is noted that transpiration also affects 
the recovery factor; however, it can be observed that the effect of 
transpiration is not too large. In fact, the transpiration effect is 
much less than 1 the apparent effect of Reynolds number for the zero 
transpiration case. Since it is known from experiment (ref. 10) that 
the recovery factor at zero transpiration is relatively independent of 
Reynolds number, doubt is shed on the validity of the recovery-factor 
determinations of this analysis. From a consideration of relative 
results, all that can be concluded is that transpiration probably does 
not affect the recovery factor by a large amount and that a value of 
r = 0.9 should suffice in design calculations. 
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In applying the results of this analysis as expressed in figures 8 
and 9, it is convenient to express the parameter £ in terms of the 
temperatures controlling a given cooling problem. If an effectiveness 
of cooling is defined by 


E 



( 12 ) 


for a surface where a heat balance (eq. (9)) applies, it is possible to 
relate C and E by 


E = 


i + 


I 

St 

(cf/2) 


( 13 ) 


Thus, once E is specified for a given design, it is possible to 
determine £ and, consequently, c f , St, and p w v w for given flight 

conditions and surface temperature. 


CONCLUDING REMARKS 


It can be concluded from the results of the present analysis that 
a transpiration cooling system is highly effective, even at high flight 
speeds. Not only does the coolant absorb a maximum amount of heat by 
attaining as its end temperature the temperature of the surface that it 
is cooling, but it has the further advantage that the effect of the 
transpiration greatly reduces both the skin friction and the amount of 
heat entering the surface. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 12, 195^ 
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APPENDIX A 


DETAILS OF ANALYSIS 


When it is assumed that changes in the x direction are negligible 
compared to changes in the y direction, equations (l) through ( 3 ) 
become 


_d_ 

ay 


p T V T 

) = 0 

(Al) 

-( 

\ £ -p—') 

(A2) 

dy N 

V dy J 



/— — if nT 3 s \ a / - cLT — du 

(pv + p'v') — c-nT + — ) = — k — + |!U — 


dy 


2 / dy V dy dy 


|iu — - pcpV'T' + upu'v' 


(A3) 


These equations can be simplified by integrating equation (Al) . The 
integral of equation (Al) is 


pv + p T v ! = constant 
At the surface of the plate 

pv = p w v w 

and 

p'v ? = 0 

Therefore, 


pv + p’v' = p w v w 

Equations (A2) and (A 3 ) thus become 

du d f - du 


p W Vw dy dy V ^ dy " PU V 


(A4) 

(A5) 

(A6) 

(A7) 


(A8) 
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PttVtt lL ( c-qT + — \ = — ( k SiZ + jlu Zii - pcpV'T’ - upu'v'^ (A9) 
W dy V P 2 / dy V dy dy y J 

It is convenient, at this point, to introduce the concepts of eddy 
viscosity, eddy thermal conductivity, and turbulent Prandtl number. The 
eddy viscosity and eddy thermal conductivity are defined as: 



(A10) 


and 


e H 


pcpV'T' 



(All) 


The usual definition of Prandtl number is pcp/k, where the fluid 
properties are based on the molecular transport of momentum and energy. 
In this analysis there is also employed the analogous quantity £m c p/ € H 

in which the properties, in effect, are based on the turbulent transport 
of momentum and energy. This quantity is called the turbulent Prandtl 
number and is given the symbol Pr^. 


When equations (A10) and (All) are introduced into equations (A8) 
and (A9) there results 


du 

Pu v w ^ - 



(l* + e M> 


du 

dy _ 


(A12) 


Pw v w 


_d_ 

dy 



d_ Z' JL + e M \ _d_ 
2 ) dy \ Pr Pr t / dy 


(c p T) + (p + e M ) 



dy 

(A13) 


Note that the bars representing time -average 
in these equations as all the terms are mean 
terms no longer appearing. 


quantities have been dropped 
values , the fluctuating 
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I.- SKIN-FRICTION RELATIONSHIP FOR Pr = Pr t = 1 
Temperature-Velocity Relationship 


For the case where Pr = Pr^ = 1, equation (A13) becomes 


PW V W 


dy 


CpT + — 
* 2 


d_ 

dy 


/ x a 

( d + €jyj) — 

dy 




(All.) 


Equation (Al4) has the identical form of equation (A12) . Thus, if u 
is a solution of equation (A12) , then 


c p T + 



au + b 


(A15) 


is a solution of (Al4) since it is linear. The terms a and b are 
constants to be evaluated at the boundary conditions: 


T = Ti at u = u x (outer edge) 


T = Ty at u = 0 (surface) 

When these boundary conditions are employed, equation (A15) becomes 


where 


— = 1 + Bu - A 2 u 2 
T W 



(A16) 


(A17) 



(A18) 


(A19) 


Because of equation (Al6) , it is only necessary to integrate the 
momentum equation (A12) to obtain both the velocity and temperature 
distribution in the boundary layer. To perform this integration, it is 
assumed that the boundary layer is divided into two parts: a sublayer. 
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adjacent to the surface, where momentum transfer is by molecular means 
alone; and an outer portion of the boundary layer, where momentum 
transfer is by eddying motion alone. 


Velocity Distribution in the Sublayer 


In the sublayer, equation (A12) is written as 


du 

P W V W ~ — 


ay 


_d_ 

ay 



(A20) 


Since p w v y is not a function of y, equation (A20) can be integrated 
directly to yield 


P 


r tVt t u = u. — + constant 
W V. * dy 


(A21) 


The constant in equation (A2l) can be evaluated by employing the boundary 
conditions 


y = o 


u = 0 


■w 


du 

4 ay 


(A22) 


Equation (A2l) becomes 


au 

V- ^ = T w + p w v w u 


Integration of equation (A23) yields 


-u 


y = 


p.du 


Jo p w v w u + T w 


(A23) 


(A24) 


It was shown in the Discussion that the extent of the sublayer under the 
conditions of blowing is not well known. In view of this, extreme rigor 
in the integration of the right member of equation (A24) is not mandatory. 
For simplicity, then, let q = My and the resulting integration yields 



PiUi 



(A25) 
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where 



(A26) 



(A27) 


ui /%r 
= \T W / 


(A28) 


Velocity Distribution in the Outer Turbulent Portion 
In the outer turbulent portion, equation (A12) is written as 


P 


W V W 


du _ _d_ 
dy dy 



(A29) 


By employing Prandtl's mixing-length concepts = pi 2 

l = Ky, equation (A29) can be rewritten in the form 



and using 


n v du d_ 

P » Tw dy - dy 



(A30) 


Because Py v y is not a function of y, equation (A 30 ) can be integrated 
directly to yield 


fty v w u = P^y 2 



+ constant 


(A31) 


At the interface of the sublayer and outer turbulent portion, it is 
required that the velocity on the laminar side match the velocity on the 
turbulent side. In addition, a matching of the laminar and turbulent 
shears is also required. By comparison of corresponding terms of equa- 
tions (A 2 l) and (A31), it is apparent that the matching of the velocities 
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and shears requires the constants to be the same. Thus, equation (A3l) 
can be rewritten as 


P W V W U + T w 



Integration of equation (A32) results in 


y = y a ex P 



Kp 1/g du 

Av + T 


W 


(A32) 


(A33) 


It will be convenient later to represent the exponential term in 
equation (A33) by the symbol g, thus 

y = y a s (A34) 

Because the primary concern of this analysis is the determination of 
skin friction and heat transfer at the surface, it is not necessary to 
perform the integration indicated in equation (A33) at this stage of the 
analysis . 


Substitution of Velocity Distribution Into 
von Karman Momentum Equation 


The von Karman momentum-integral equation for the case with surface 
blowing at constant pressure can be written as 


T w + p w v w u i = ^ 


pu (ux-u) dy 


(A35) 


Equation (A35) can be made dimensionless by dividing through by Piui 2 . 
It is also convenient to replace the variable of integration y by u. 
The resulting equation is 


+ Vw = _L 

2 PiUi dx 



(A36) 


The integral in the term on the right represents the momentum thickness 9 . 
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To evaluate the momentum thickness rigorously in equation (A 36 ), it 
would be necessary to divide it into two parts. The first part, with 
the lower bound as 0 and the upper bound as u a , would represent the 
contribution of the sublayer. The second part, with the lower bound as 
u a and the upper bound as 1 , would represent the contribution of the 
outer turbulent portion. In the first part (dy/du) would be evaluated 
from equation (A 23 ), while in the second part it would be evaluated from 
equation (A33)- For the incompressible case without blowing, these 
integrations can be performed very simply in closed form. It is found 
that this two-part system can be approximated very closely (less than 
1-percent difference in Cf(Rg) curves) by the integral of the second part 
in which the lower bound is set equal to zero. It would not be expected 
that the introduction of compressibility and blowing should markedly 
alter this behavior; consequently, the simpler one-part system will be 
used in the present analysis to determine 6 . 

When (dy/dfr) is evaluated from equation (A33)> the expression for 
momentum thickness becomes 


e = Ky a 


£X*r 


u 


(l-uV 


J 


du 


Pi ui 




u + 


noting that 



(A37) 


(A38) 


Equation (A37) cannot be solved in closed form; however, if the condition 
of compressibility is removed, a closed form solution for Q can be 
obtained easily. It is found from the solution for the incompressible 
case that the first term resulting from an integration by parts yields 
results which have the same form as the results of the complete integra- 
tion. Since we are primarily concerned with the form of the solution, 
absolute magnitudes can be adjusted by slight changes in u a and K which 
are arbitrary constants ; only the first term resulting from an integration 
by parts of the compressible case will be used in this analysis. 


To perform the integration of equation (A37)> 
replace u as the variable of integration by g. 

-J s / ru 

dg = — — g du 


K 


V U + 2 


it is necessary to 
From equation (A 38 ) 


(A39) 
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Substitution of equation (A39) into equation (A37) yields 

g(i) 


9 = y a / 

g(o) 


3 Q u (l-u) dg 


(a4o) 


Because the pressure across the boundary is assumed to be constant ^ the 
perfect gas law and equation (Al 6 ) indicate that 




(A4l) 


Equation (A40) then becomes 


(©/ 


g«. 


u(l-u) dg 

9 ■ yal:f W^4(o) 


(AU2) 


When equation (A42) is integrated by parts, there results 


0 = y a =r- 


Ti 


■W- 


u=i 


u(l-u)g 


1+Bu-A 2 ^ „ 


u=o 


_r 1 l-2u+(A a -B)u g du 

J 0 8 ( 1+Bu-A 2 iI 2 ) 2 


(A43) 


From an examination of equation (A38) it is apparent that both g(l) 
and g(0) are finite; consequently, the first term on the right of 
equation (A43) vanishes at the bounds. If g is again substituted for 
ft as the variable of integration, equation (A43) becomes 


e = 


r 

K 


g(°) 


. , / P W V W ~ 

S( } {l-2u+(A 2 -B)u 2 l«/ P X U! U + 2 

( l+Bu-A 2 ^ 2 ) 3/2 


dg (A44) 


Again, integrating by parts results in 


0 = - 


K 


1 


-1 / P W V W ~ ff 
)u^l V PiUi 2 


U=1 


f-2u+(A 2 -B)u 2 | P 1 U 1 

(l+Btf-A 2 ^) 372 


u=o 


(a4 5 ) 


+ additional terms 
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When the additional terms are omitted and the limits of integration are 
introduced 


e = 


ya ^) / 

K 


g(i) J¥ 


p w v w 

Pl u l 


(i+b-a 2 ) 1/s 

For practical cases g(l) » g(0); therefore, 

X 1 / 2 


Ya 


0 = 


1 

T W- 

K 


J¥ 




(A k6) 


P W V W 

Pxu x 


(l+B-A 2 ) 1 / 2 


g(l) 


(A47) 


But 

T W 


(l+B-A 2 ), so that 


0 




g(l) 


Employing equation (A26) 



g(D 


Equation (A49) can be manipulated into the form 


(A48) 


(AA9) 



(A50) 


Before integrating the right member of (A50), it is preferable to obtain 
the relationship of Cf in terms of the length Reynolds number. 
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Equation (A36) can be rewritten as 



(A51) 


when equation (A^9) is employed. On integration, neglecting the effect 
of laminar sublayer 


R x = 


C -P 

,-r < x > 


Cf ^ 

> 00 


I A- 


( p w v w/ PlUl) 

( c f/2) 


g(D 


_£ Vw 

- + p x u x 


(A52) 


Integration by parts yields 


Rx = 


’i./‘ - ■»' 

(c f /2) + ( p w v w/ p i u i) 


(x) 


+ additional terms 


(A53) 


f ^ 


For large values of R x the first term predominates. Also g(l) is 
finite ; therefore , 


r x = ^ 

x K 


g(l) 


_£ (-L + PwV w 
2 V 2 p iU;l 


or 


VTW^ 


PttV, 




Rv . 


pV/2 


K' « 


/ 


P W V W ~ ff 

P1U1 + 2 


(A54) 


(A55) 


It is apparent that the right members of equations (A50) and (A55) are 
identical. This does not mean that the left members are equal, but just 
that the K and u+ (also y+) used in each of the equations are different 
(see Discussion) 
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The integral appearing in equations (A50) and (A55) can be written 


as 


K 0r' 


iY /a r 1 
V 


/ 


du 


a 


V PjUj. 2 


(A56) 


when equation (Al6) is employed. The integral (A56) can be rewritten as 


where 


K 


AiY /2 

VT„/ J x 


du 


/ Y- V .y O' (u+ct ) (u-a) ( P-u) 

y Pi u i a 


a = 


P = 


B - J B 2 +4A 2 
2A 2 

B + B g +4A 2 
2A 2 


a = 


(c f /2) 


P W V W 

Pl u i 


(A57) 


(A58) 

(A59) 

(A60) 


Letting 


sin 2 9 = c (p-u) 


(A6l) 


then du = - — sin cp cos cp dcp and the integral (A57), on transforming 
the variable of integration from u to cp, becomes 

\l/2 


2Kl .|) p*(D 


/ 


— cos cp dcp 

A L 


J Vine) J (p+a) - i sin 2 cp (P-a) - i sin 2 cp 


(A62) 


If c is set equal to then integral (A62) becomes 
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2K 


M 1/2 

V 

P w v w ^ 


/VV 

V P1U1 



dcp 


3+o 
3 -a 


. P 

sin 


9 


Since Legendre’s standard form of the elliptic integral of the 
is 


F(k, q>) = 


^9 


dcp 


J 1-k 2 


sin 2 cp 


Integral (A 63 ) can be rewritten as 


2K [ iV 2 
V 


/vw^ 

y p x u x 


F [km, 9 m (u a )] -F[km, q>(lj] 


where 


and 


v2 - p±2 ^ l 

^ ■ 3-a < 


sin 2 <P m (u) 


3-u 

3+c 


If c in equation (A62) is set equal to l/3-a, integral (A63) 


/tA 1 / 2 
2K 1 tV 


j W 

V Px u i 


V3 +cr 


{ F [ k n> 9 n (u a )] - F [k n , <p n (l)] 


where 


k 2 - iz® 

^ ■ 3^ K1 


• 2 ^ / ~ \ 3-U 

sin ^ (u) - 


(A63) 


first kind 

(A64) 


(A 65 ) 


(A66) 


becomes 

(A67) 


(a68) 
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II.- SKIN FRICTION - HEAT TRANSFER RELATIONSHIP FOR Pr / 1. 


This part of the analysis begins with the integration of equations 
(A12) and (A13), resulting in 


p w v w u = < p + € m) C 1 


(A 69 ) 


and 


P W V W l C P T + 2 


^ + ( p+e M) 


Pr Pr^ 


dCr,T 


<*y 


s£) 

dy 


+ Co (A70) 


where Ci and c 2 are arbitrary constants. At the surface y = 0, we 
have 


U = 0 
T = T W 


du 

T w=^ 


, dT 

H W = - k — 


e M = 0 


W 

IT 

dy 


W 


jx_ dc p T 
Pr dy 


W 


At the surface, equations (A 69 ) and (A 70 ) thus become 


0 = T w + Cl 


0w v w c P% = - Iw + c 2 


(A71) 

(A72) 


When the arbitrary constants evaluated from equations (A 71 ) and (A72) are 
substituted into equations (A 69 ) and (A 70 ), there results 


du 


P W V W u + t w = ( P + e M ) ^ 


'W 


dy 


(A73) 


p w v w 


( T _ T ) . Hi 1 _ q = fjL + + (n+e M )u — (A7U) 

)U V + 2 4 w \Pr Pr t / dy ^ M dy 
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Dividing the terms of equation (A74) by the terms of equation (A73) 
yields 


P W V W 


[c p (T-T w ) + Jjf] - q* 

P w v w u + T w 


( JL . e M ^ 

V Pr + Pr t y 1 dc p T < ^ 

( U + eM ) du 


(A75) 


It is convenient in the processes which follow to introduce the following 
symbols : 

F = ^ 

Pl u i 


E = c p (T-T w ) + ^ 

Equation (A75) can then be written as 


P = 


( P + e M ) 


M_ fc M 
Pr + Pr + 


S = 


% 

P,u, 


dE 

du 


PF 

Fu + Cj./2 


E = (l-P) U! 2 U - 


PS 

Fu + c-f/2 


In the range 


0 <u <u a P = Pr 

u a <u<l P = Pr t 


(A76) 


Equation (A76) must then be integrated in two stages with the results 
combined at the interface where ft = u a . The integration, with the use 
of an integrating factor, is straightforward. After some algebraic 
manipulation, there results back in the original notation 


J*w_ 

PjL U l 


Ul 


- c p Ti - r — 


Vw 

PX U 1 


C J. 

2 


\ Px u i 2 / 


-Pr /P„v, 


p w v w 

P x u i 


Pr t 

+-t ) - 1 


(A77) 


If we define the heat-transfer coefficient as 


h 


L W 


% - Tx - r 


Ul 
2 Ct 


(A78) 


the final relationship relating heat transfer and skin friction is 


St 


'f 

2 


f P W V W 
V Pi u i. 


_ / u f 


1-Pr 


( P W V W ~ 
\ Pl u l Ug 


N Pr-Pr t 
Cf ' 1 


p w v w 

p l u i 


♦sr 


- _£ 


(A79) 


The expression defining the recovery factor is 


r = 1 


p W v W ^ 

p i u iy 


+ ^ ) (2-Pr)(2-Pr t ) - (2-Pr) 


P W V W c f 
p l u i T 


PWVW S C£_ 
Pl% a 2 


2 -Pr t 


( p W v W | c f 

\PlUi 2 


Pr t 


(Pr t -Pr) 


@g) <* - PrtX* - Pr) 


(2 


- 2 


-**>(?)•* 


7p»» w~ cA Pr - prt 


7^*7 — 7 

\oTuT/ (2 • Pr ‘ ,(2 


p w v w 

PlUl 


v) 


Prt 




(A8o) 


UJ 

VJ 1 


NACA TW 33^1 


36 


NACA TN 33^1 


REFERENCES 


1. Eckert, E. R. G., and Livingood, John N. B.: Comparison of Effec- 

tiveness of Convection-, Transpiration-, and Film-Cooling Methods 
With Air as Coolant. NACA TN 3010, 1953- 

2. Eckert, E. R. G.: Heat Transfer and Temperature Profiles in Laminar 

Boundary Layers on a Sweat-Cooled Wall. AAF, Air Materiel Command, 
Wright Field. TR 5646, 1947. 

3. Brown, W. Byron, and Donoughe, Patrick L.: Tables of Exact Laminar- 

Boundary-Layer Solutions when the Wall is Porous and Fluid Proper- 
ties are Variable. NACA TN 2479? 1951. 

4. Rannie, W. D.: A Simplified Theory of Porous Wall Cooling. 

(Progress rept. 4-50. ORDCIT proj. Contract W-04-200-0RD-455? 
Ordnance Dept. Power Plant Lab. proj. MX801, Contract W-535-AC- 
20260, Air Mat. Com.) C»I.T., Pasadena (Guggenheim Aero. Lab.), 

Jet Propulsion Lab. 

5. Friedman, Joseph: A Theoretical and Experimental Investigation of 

Rocket-Motor Sweat Cooling. Jour. American Rocket Soc., No. 79? 
December 1949? pp. 147-154. 

6. Crocco, L. : An Approximate Theory of Porous, Sweat, or Film Cooling 

With Reactive Fluids. Jour. American Rocket Soc., vol. 22, no. 6, 
Nov. -Dec. 1952, pp. 331-338. 

7. Duwez, P., and Wheeler, H. L.: Experimental Study of Cooling by 

Injection of a Fluid Through a Porous Material. Jour. Aero. Sci., 
vol. 15 , no. 9? Sept. 1948, pp. 509-521. 

8. Jakob, M. , and Fieldhouse, I. B.: Cooling by Forcing a Fluid Through 

a Porous Plate in Contact With a Hot Gas Stream. Heat Transfer 
and Fluid Mech. Inst. (Paper presented at June 22-24, 1949 meeting) 

A. S . M. E . 

9. Mickley, H. S., Ross, R. C., Squyers, A. L., and Stewart, W. E.: 

Heat, Mass, and Momentum Transfer for Flow Over a Flat Plat With 
Blowing or Suction. M.I.T., Sept. 1952, NACA Contract NAw 5980. 

10. Rubesin, Morris W.: A Modified Reynolds Analogy for the Compressible 

Turbulent Boundary Layer on a Flat Plate. NACA TN 2917? 1953* 

11. Li, Ting-Yi, and Nagamatsu, H. T.: Effect of Density Fluctuations 

on the Turbulent Skin Friction on a Flat Plate at High Supersonic 
Speeds. Proc. of the Second Midwestern Conference on Fluid 
Mechanics, 1952. Ohio State Univ. Engineering Experiment Station 
Bulletin No. 149, pp. 263-275. 


NACA TN 3341 


37 


12. Coles, Donald: Measurements in the Boundary Layer on a Smooth Flat 

Plate in Supersonic Flow. C.I.T. Thesis, Pasadena, California, 

1953. 

13. McAdams, William Henry: Heat Transmission. McGraw-Hill, New York, 

I9te, P. 175. 

14. Fage, A., and Falkner, V. M. : Relation Between Heat Transfer and 

Surface Friction for Laminar Flow. R. & M. No. 1408, British 
A.R.C . , 1931. 

15* Colburn, Allan Philip: A Method of Correlating Forced Convection 

Heat Transfer Data and a Comparison With Fluid Friction. Trans. 
Am. Inst. Chem. Engr., vol. 29 , 1933, PP* 17^-210. 

16. Rubesin, Morris W. : The Effect of an Arbitrary Surface -Temperature 

Variation Along a Flat Plate on the Convective Heat Transfer in 
an Incompressible Turbulent Boundary Layer. NACA TN 23^-5, April 
1951. 


38 


NACA TN 3341 























Analytical results with" 
parameters K and u£ 
based on 
velocity survey 


NACA 




'SpK 

gf 



PI 

Mm 




• -■ • - 


Ijjpirii 




-TtTtfll! 1 

|§j|§| 






— 

















M 











In iff 



Liar 

Sjjm 


np? 





' 


gpH H$H 



If ilT 1 mTm 

r 



• • : 




-III Hill 


i! 




••f" i 


Sjj 

St ■ f -f-H rf frr 

r 

w|8 




■ 1 

11 
jmj tm 


l^rfeppl 


Hilil 



ifUl 





rfft 

!!■■ 









■ 








o Schu 

lltz- 

Grunow ! 













;;;;;;;;; 

• ITI nknui/>m 

C rr\ D r\£ rw r\ir\ r* IO U 


> H 

(em 

fffi 

pf 

is 


miuii 

gjfffl 

IfSll 


imiiiifiiiiriiiiii 



iacsip 

jE 

jgg|| 

... 



1 





ill 


:::: 







warn 



44 


It 

l H 1(11 111 










mm 







±± 





jagfjg p 


































r- 

tHt 






.... L 















JL. . 





— — 



11 

— - 






1MB 








TtftttW 








Mill 



— — J— 




IM 


mllllll 



=1 

1 

iiiii 

B li 



■ 11 



Hi-miiiii 




Figure 1.- Direct measurements of local skin friction in low-speed flow compared with analytical 
results based on alternative choices of arbitrary parameters. 
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Figure 2.- Data of reference 9 for incompressible flow, no transpiration case compared with 

results of present analysis. 
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Figure 4.- Data of reference 9 for incompressible flow with uniform transpiration 

compared with present analysis. 
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©Run 19a- Constant blowing gas temperature | 
□ Run 19b — Constant surface temperature 
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Figure Data of reference 9 for incompressible flow with uniform transpiration — = 0.006 


compared with present analysis. 
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Figure 6.- Data of reference 9 for incompressible flow with varying transpiration rate 



compared with present analysis. 
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Figure 8. - 


(a) M = 0; Ty = free-stream stagnation temperature. 

Effect of transpiration rate on local skin-friction coefficient. 
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Figure 8.- Continued 
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(c) M = 2; Ty = free-stream stagnation temperature. 


Figure 8.- Continued. 
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Figure 8.- Concluded. 
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(a) M = 0; T w = free-stream static temperature. 

Figure 9 *~ Effect of transpiration rate on local ratio of Stanton number to skin-friction 
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(b) M = 2; Ty = free-stream static temperature. 
Figure 9* - Continued. 
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(c) M = 2; Ty = free-stream stagnation temperature. 
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Figure 9*- Continued 
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Figure 9 .- Concluded 
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